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2  Velocities  caleulated  from  MULTIA:  (a)  after  6000  s  and 
(b)  after  4000  s  The  extent  of  the  mushy  zone  is  also 
shown.  Dimensions  are  in  m. 
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is  800  pm  and  (b)  primary  arm  spaeing  is  300  pm.  The  extent 
of  the  mushy  zone  is  also  shown.  Dimensions  are  in  m. 

4  Maerosegregation  of  oxygen  caleulated  from  MULTIA 

after  8700  s  of  ingot  solidifieation  with  a  primary  arm  spaeing 
of  300  pm:  (a)  element  size  is  14.3  mm  x  21.7  mm  and 
(b)  7.2  mm  x  10.8  mm.  The  extent  of  the  mushy  zone  is  also 
shown.  Dimensions  are  in  m. 

5  Maerosegregation  of  oxygen  caleulated  from  MULTIA 
after  16000  s  of  solidifieation  with  a  primary  arm  spaeing 
of  300  pm:  (a)  element  size  in  mesh  is  14.3  mm  x  21.7  mm 
and  (b)  7.2  mm  x  10.8  mm.  Dimensions  are  in  m. 

6  Velocities  caleulated  from  MULTIA  with  a  primary 

arm  spaeing  of  300  pm  and  element  size  of  7.2  mm  x  10.8 
mm:  (a)  after  2000  s  and  (b)  after  4000  s  of  ingot  solidifieation. 
The  extent  of  the  mushy  zone  is  also  shown. 

Dimensions  are  in  m. 

7  Velocities  caleulated  from  MULTIA  with  a  primary 
arm  spaeing  of  300  pm  and  elements  of  7.2  mm  x  10.8 
mm:  (a)  after  6000  s  and  (b)  after  8700  s.  The  mushy 
zone  is  also  shown.  Dimensions  are  in  m. 


Summary 

Simulations  of  VAR-ingots  of  Ti-6A1-4V,  with  emphasis  on  predicting 
macrosegregation  of  oxygen,  are  presented.  Processing  conditions  for  the  simulations  were 
selected  in  with  RMI  Titanium  of  Niles,  Ohio,  who  provided  the  output  from  their  software 
(BAR),  which  was  used  to  assist  in  analyzing  an  experimental  ingot.  BAR  does  not  capture 
macrosegregation,  so  its  calculated  Lorentz  forces  were  extracted  and  included  as  body  forces 
along  with  the  gravity  body  forces  in  MULTIA.  MULTIA  simulations  showed  that 
segregation  patterns  in  the  ingots  result  largely  from  the  strong  convection  in  the  melt  pool. 
Permeabilities  in  the  mushy  zones  were  less  important.  Simulations  revealed  more  details  in 
the  macrosegregation  patterns  when  the  element-size  selected  for  the  calculation  was  reduced 
to  7.2  mm  x  10.8  mm.  The  extent  of  the  segregation,  as  the  difference  between  the  maximum 
0.22  wt.%)  and  minimum  (=;  0.15  wt.%)  concentrations  of  oxygen  were  about  the  same  for 
both  element  sizes.  A  band  of  negative  segregation  near  and  running  parallel  to  the  ingot- 
surface  was  simulated,  but  the  finer  sized  elements  revealed  a  more  intense  band.  Another 
feature  was  a  broad  vertical  band  of  positive  segregation,  which  originated  at  the  center- 
bottom  and  migrated  to  the  half- width  position,  where  it  persisted  to  about  5/8  of  the  ingot 
length.  Since  primary  arm  spacings  in  titanium  alloys  are  not  available,  primary  dendrite  arm 
spacings  in  Ti-6A1-4V  were  estimated.  A  summary-status  of  the  use  of  software  by  VAR 
titanium- ingot  producers  in  the  USA  is  also  given.  In  its  initial  stages,  the  research  was  on 
"Formation  of  Microporosity  in  Nickel-Base  Thin  Wall  Castings,"  so  a  brief  summary  of 
microporosity  in  experimental  castings  is  presented  (Appendix  C). 

I.  Introduction 

In  order  to  simulate  macrosegregation  in  VAR  ingots,  an  “in-house”  code,  MULTIA, 
was  used.  It  is  based  on  a  complete  model  of  transport  phenomena  (energy  equation, 
conservation  equations  for  each  solute,  momentum  equation,  and  continuity).  The  mushy  zone 
is  treated  as  a  porous  medium  in  which  the  permeability  varies  with  fraction  solid. 
Mathematical  details  underlying  MULTIA  are  in  Appendix  A.  Convection  in  the  melt  pool  is 
driven  by  both  buoyancy  and  Lorentz  forces.  The  latter  are  particularly  important  in  Ti-ingot 
processing.  Selection  of  the  processing  conditions  for  the  simulations  was  done  with  RMI 
Titanium  of  Niles,  Ohio,  who  provided  the  output  from  a  proprietary  code,  BAR.  BAR 
simulates  convection  in  the  melt  pool  and  demarcates  the  mushy  zone.  A  first-order 
turbulence  model  with  Lorentz  forces  is  included,  but  BAR  is  not  set  up  to  capture 
macrosegregation.  In  this  work,  Lorentz  forces  were  extracted  from  BAR  and  included  as 
body  forces  in  MULTIA  so  that  the  interaction  between  the  melt  pool  and  the  mushy  zone 
were  simulated.  Equilibrium  partition  ratios  of  alloy  elements  were  provided  by  the  Air  Force 
Materials  Laboratory  (AFML). 

In  this  report,  simulations  of  the  macrosegregation  in  VAR-ingots  of  Ti-6A1-4V  are 
reported.  The  output  comprises:  the  temperature  field  and  the  position  of  the  mushy-zone;  the 
concentrations  of  Al,  V,  O  and  Fe  and  their  variations  across  the  diameter  of  the  ingot;  and 
velocities  in  the  melt  pool  and  the  mushy  zone.  Subsidiary  information  includes  the  Lorentz 
forces  in  the  melt  pool,  the  volume  fraction  of  solid  throughout  the  mushy-zone,  and 
microsegregation  characteristics. 
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II.  Estimating  Dendrite  Spacings  in  Ti-6A1-4V 

A.  Connection  between  Permeability  and  Dendritic  Structure. 

Since  permeability  has  units  of  ,  a  length  scale  of  the  dendritic  structure  must  be 
assigned  to  estimate  the  permeability  and  hence  the  strength  of  convection  during 
solidification  of  alloys.  We  devise  relationships  for  estimating  the  dendritic  spacings  in  Ti- 
6A1-4V  alloy  for  the  purpose  of  simulating  the  formation  of  macrosegregation  in  ingots  of  the 
subject  alloy.  The  convection  in  the  mushy  zone  is  modeled  with  a  momentum  equation, 
which  includes  inertial  and  viscous  terms  along  with  the  Darcy  term  [1],  in  which  the 
permeability  tensor  appears.  Permeability  ( K )  has  dimensions  of  length-squared.  For 
columnar-dendritic  growth  [2]: 


K  =  ^,^f,(g,)  (1) 

where  fjCg^)  is  a  function  of  the  volume  fraction  of  the  interdendritic  liquid  and  is  the 

primary  dendrite-arm  spacing.  Dendrite  spacings  are  reported  for  aluminum  alloys,  ferrous 
alloys  and  Ni-base  alloys  [4,  5],  but  relatively  few  data  on  titanium  alloys  exist.  Nurminen 
and  Brody  [6]  measured  and  reported  secondary  arm  spacings  ( d2,  pm )  in  several  binary 
alloys,  including  Ti-Al  and  Ti-V.  Data  for  primary  arm  spacings,  however,  are  not  available. 
Here  the  secondary  arm  spacings  for  Ti-Al  and  Ti-V  alloys  are  summarized,  from  which 
primary  arm  spacings  are  estimated. 

B.  Secondary  Dendrite  Arm  Spacings 

Nurminen  and  Brody  [6]  measured  and  correlated  secondary  dendrite  arm  spacings 
with  solidification  time  (^,  s )  as 


d^  =  CO"  (2) 

where  d2  is  the  secondary  spacing  (pm)  and  C  and  n  are  constants.  The  constants  for  the  Ti-V 
and  for  the  Ti-Al  alloys  are  about  the  same,  so  each  was  treated  as  one  data  set.  C  is 

C  =  184.36  -  41.453  X  +  4.2375  x"  -  0.1577  x'  (3) 

where  x  is  either  wt.%  A1  or  wt.%  V.  Ripening  theory  applied  to  secondary  dendrite  arm 

spacing  predicts  n  of  1/3  [4].  The  values  found  range  from  0.312  to  0.444,  in  rough 

agreement  with  the  theoretical  value.  As  a  function  of  concentration, 

n  =  0.3411  +  0.0053  X  (4) 

Since  the  binaries  have  similar  constants,  Eqs.  (3)  and  (4)  were  utilized  by  taking  x  as 
the  sum  of  the  concentrations  of  A1  and  Ti  in  the  commercial  alloy,  Ti-6A1-4V.  For  the  Ti64 
alloy,  we  use  a  composition  (in  wt.%)  of  Ti  -  5.72  A1  -  4.10  V  -  0.22  Fe  -  0.183  O;  hence  C  = 
36.585  and  n  =  0.3931  from  Eqs.  (3)  and  (4),  respectively. 

In  a  rapidly  solidified  Ti-alloy  (Ti-15V-3Al-5Sn-3Cr),  a  relationship  between  the 
secondary  dendrite  arm  spacing  and  cooling  rate  is  presented  [7].  This  relationship,  however. 
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is  based  on  applying  a  heat  transfer  model  for  a  rapid  solidification  scenario;  the  empirical 
results  of  Nurminen  and  Brody  [6]  are  deemed  more  reliable,  especially  for  cooling  rates 
encountered  in  ingots  and  castings. 

C.  Primary  Dendrite  Arm  Spacings 

In  order  to  estimate  primary  dendrite  arm  spacing  in  Ti64,  the  data  collected  by 
Whitesell  [5]  for  both  primary  and  secondary  dendrite  arm  spacings  in  Ni-base  superalloys 
were  consulted.  In  terms  of  a  cooling  rate,  the  dendrite  arm  spacings  are  approximated  by 

d,  =  160  (5) 

and 

=  160^,"“^‘'‘  (6) 


The  cooling  rate  is  in  K/s,  and  the  arm  spacings  are  in  pm.  The  cooling  rate  in  Eqs.  (5)  and  (6) 
are  St  =  GR,  where  G  and  R  are  the  thermal  gradient  (K/cm)  and  solidification  rate  (cm/s) 
measured  in  the  liquid  at  the  dendrite  tips.  Since  this  G  differs  from  that  through  the  mushy 
zone,  then  St  differs  from  the  usual  cooling  rate  based  on  solidification  time,  namely: 


where  Tl  is  the  liquidus  and  Te  is  the  non-equilibrium  solidus. 

A  relationship  between  St  and  s  was  sought  by  making  use  of  empirical  cooling 
curves  in  a  Ti-alloy  casting  given  by  Nurminen  and  Brody  [6]  and  cooling  curves  for  Fe  -  26 
wt.%  Ni  given  by  Flemings  et  al.  [8].  The  ferrous  alloy  has  a  liquidus  lower  than  that  of  the 
Ti64  alloy  (1653  °C),  but  it  is  similar  to  the  Ti64  alloy  in  that  it  has  a  relatively  small 
solidification  temperature  range.  In  all  eight  sets  of  cooling  rates  were  regressed,  resulting  in 

f  =  6.614  (8) 


The  cooling  rates  covered  by  Eq.  (8)  are  up  to  5.3  K/s  for  s  and  0.75  K/s  for  St. 

To  convert  from  cooling  rate  to  solidification  time,  solidification  temperatures  of  the  Ti64 
alloy  of  the  composition  Ti-5.72  Al-4.10  V-0.22  Fe-0.183  O  were  used.  The  liquidus  is  1926 
K.  The  temperature  at  the  end  of  solidification  must  lie  between  the  equilibrium  solidus 
(1964.5  K)  and  the  non-equilibrium  solidus  assuming  no  diffusion  in  the  solid  during 
solidification  (1881.4  K).  Using  Eq.  (2)  with  the  constants  for  the  Ti64  alloy,  solidification 
times  were  varied  from  1  to  100  s,  and  two  cooling  rates  in  the  mushy  zone  were  calculated 
for  each  solidification  time.  For  one  cooling  rate  it  was  assumed  that  Tg  was  1964.5  K,  and 
for  the  other  it  was  1 88 1 .4  K.  Then  the  cooling  rates  were  used  to  deduce  the  primary 
dendrite  arm  spacings  by  the  ratio  given  by  combining  Eqs.  (5)  and  (6): 
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(9) 


4.893 


-0.03 


Finally  the  dendritic  spacings  were  put  in  terms  of  s,  which  is  the  usual  cooling  rate  given  by 
Eq.  (7). 

Dendrite  spacings  calculated  by  the  above  procedures  are  shown  in  Fig.  1 .  The 
secondary  dendrite  arm  spacing  is  the  same  as  Eq.  (2)  with  C  =  36.585  and  n  =  0.3931.  There 
are  two  lines  (practically  superimposed)  for  the  primary  dendrite  arm  spacing,  depending  on 
the  selection  of  the  temperature  for  the  solidus.  The  average  of  the  two  lines  is 

d,  =  168.62  (10) 


Fig.  1 .  Calculated  primary  and  secondary  arm  spacings  in  Ti-6A1-4V. 
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III.  Macrosegregation  Simulations 

A.  Role  of  Primary  Dendrite  Arm  Spacing  on  Macrosegregation 

Simulations  of  solidification  of  an  ingot  of  Ti-6A1-4V  with  a  melting  condition 
provided  by  RMI  Titanium  Company  (Proposed  Case  11).  Two  ingots  are  simulated;  one 
is  simulated  assuming  a  primary  arm  spacing  of  800  pm  and  the  other  with  a  primary  arm 
spaeing  of  300  pm.  By  studying  these  dendritic  length  seales,  the  role  of  permeability  on 
the  overall  maerosegregation  was  ascertained.  It  is  shown  below,  that  with  the  larger 
dendrite  arm  spacing,  the  simulated  maerosegregation  is  somewhat  stronger,  but  the 
differenee  is  not  as  great  as  expeeted  if  flow  throughout  the  mushy  zone  was  important. 
Therefore,  most  of  the  maerosegregation  is  the  result  of  the  interacting  flows  between  the 
liquid  pool  and  the  leading  part  of  the  mushy  zone. 

B.  MULTIA  Simulation  of  the  VAR-Ingot 

In  VAR-processing,  the  Lorentz  force,  the  gravity  foree  {i.e.,  buoyancy  force)  and 
the  solidification  shrinkage  drive  the  eonvection.  The  pool  flow  of  the  VAR  ingot  is 
eaused  by  a  eompetition  between  the  Lorentz  force  and  buoyancy  force.  The  buoyancy 
force  is  dominant  for  low-current  VAR-ingots  {e.g.,  Ni-base  superalloys  [9]),  whereas  the 
Lorenz  foree  is  dominant  for  the  high-eurrent  VAR  ingot  (e.g.,  Ti-base  alloys). 

The  MULTIA  model  ineorporates  solidifieation  theory,  mainly  to  aceount  for  the 
partitioning  of  each  alloy  element  between  solid  and  liquid.  The  density  of  the  liquid  in 
the  mushy  zone  or  in  the  all-liquid  zone  depends  on  its  constitution  and  temperature; 
hence  the  liquid  density  varies  spatially  and  temporally  and  is  subjeet  to  gravity-driven 
eonvection  [10-13]. 

The  program  starts  with  a  mesh  that  has  four  layers  of  elements  in  the  vertieal 
direetion.  As  mass  enters  the  liquid  pool  from  the  electrode,  new  layers  of  elements  are 
added  at  the  top  of  the  ingot.  The  addition  of  layers  stops  when  the  maximum  height 
speeified  is  reaehed.  The  dimension  of  the  fully-grown  domain  is  43  em  in  width  and 
165  em  in  height.**  There  are  30  elements  across  the  width  and  76  elements  along  the 
height  of  the  complete  ingot,  and  eaeh  element  in  the  MULTIA-mesh  is  14.3  mm  x  21.7 
mm. 

As  presently  struetured,  MULTIA  is  not  furnished  to  ealculate  the  thermal 
boundary  conditions  and  electromagnetic  forces  in  the  VAR  proeess.  Therefore,  to  effeet 
a  MULTIA-simulation,  the  thermal  boundary  conditions  and  the  eleetromagnetic  body- 
forces  (Lorentz  forees)  were  taken  from  the  BAR-output.  For  convenience  and  to  avoid 
interpolations,  the  mesh  of  the  MULTIA  simulation  is  constructed  to  be  identical  to  that 
of  the  BAR  simulation.  From  the  BAR  results,  the  heat  transfer  coeffieients  on  each 
boundary  were  obtained  and  applied  as  thermal  boundary  eonditions  for  a  MULTIA 
simulation.  At  the  boundary  nodes  along  the  left  side  (line  of  symmetry),  no  thermal  flux 
is  employed.  The  MULTIA  simulation  starts  with  an  all-liquid  alloy  of  the  nominal 
eomposition  at  a  uniform  temperature  of  1950  K  (17  K  higher  than  the  liquidus)  in  the 
bottom  four  layers  of  elements. 

Figure  2  shows  the  results  from  the  MULTIA  simulation  and  illustrate  the 


Cartesian  coordinates  are  used  for  the  results  reported  here.  Axi-symmetric  coordinates  can  be  used,  but 
the  computer  runs  are  longer  than  when  Cartesian  coordinates  are  used. 
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calculated  velocities  in  the  liquid  pool.  These  figures  represent  output  after  6000  s  and 
8700  s  of  solidification  of  the  ingot.  They  show  conveetion  patterns,  whieh  are 
qualitatively  similar  to  the  BAR  result.  The  maximum  velocity  in  the  MULTIA-result 
(~40  cm/s)  is  greater  than  that  in  the  BAR-result  (~17  em/s).  In  MULTIA  the  viseosity 
of  the  liquid  equals  the  “moleeular”  viseosity,  whereas  BAR  employs  an  “eddy”  viseosity 
to  aeeount  for  turbulence. 

Figures  3(a)  and  3(b)  show  the  resulting  macrosegregation  of  oxygen  after  8700  s 
of  solidifieation  of  the  ingot.  Figure  3(a)  is  the  result  simulated  with  a  large  primary  arm 
spaeing  of  800  pm  for  the  permeability  in  mush  zone,  and  Fig.  3(b)  is  with  a  small 
primary  dendrite  arm  spaeing  of  300  pm.  Notice  that  there  is  a  band  of  positive 
segregation  of  oxygen  near  and  at  the  mid-width  of  the  ingot  for  the  both  eases.  The 
simulation  with  the  larger  primary  spaeing  enhances  segregation.  The  0-concentrated 
region  in  the  upper  area  near  the  wall  is  redueed  in  the  ease  with  the  smaller  arm  spacing. 
For  the  larger  spacing  (800  pm),  the  maximum  eoncentration  of  oxygen  is  0.334  wt.%, 
which  is  83  %  higher  than  its  nominal  eomposition  (0.183  wt.%  O).  For  the  smaller 
primary  arm  spaeing  (300  pm),  the  maximum  eoncentration  of  oxygen  is  0.280  wt.%, 
which  is  53  %  higher  than  its  nominal  composition. 

Near  the  right  side  of  eaeh  simulated  ingot,  there  is  a  vertieal  band  depleted  in 
oxygen.  Oxygen  has  a  partition  ratio  greater  than  one,  whereas  iron  has  a  partition  ratio 
less  than  one.  The  bands,  therefore,  are  enriched  in  iron  and  depleted  in  oxygen;  the 
bands  appeared  to  be  weak  freekles.  As  diseussed  in  the  next  section,  however,  the 
freekling-mechanism,  in  which  less  dense  liquid  in  the  upper  part  of  the  mushy  zone  near 
the  dendrite  tips  flows  vertieally  upward  and  exits  into  the  liquid  pool,  was  discounted. 

C.  Sensitivity  of  Calculations  on  Mesh  Spacing 

The  sensitivity  of  finite  element  mesh  spaeing  on  MULTIA-simulations  of 
maerosegregation  was  studied.  A  systematie  sensitivity  analysis  of  mesh  spaeing  on 
simulating  maerosegregation  during  directional  solidification  (DS)  [14]  was  consulted. 
The  recommended  elements  should  be  no  larger  than  2  d\  in  the  horizontal  direction  and 
1.5  DIV  in  the  vertical  direction),  where  d\  is  the  primary  arm  spacing,  D  is  the 
diffusivity  in  the  liquid,  and  V  is  the  growth  rate.  A  sensitivity  analysis  of  mesh  spacing 
was  applied  to  a  DS  domain  of  Ti-6A1-4V  alloy  with  goal  of  finding  an  optimum  mesh 
that  could  be  applied  to  simulations  of  VAR-ingots.  No  freekles,  however,  were  found  in 
the  Ti-6A1-4V  alloy,  and  in  this  DS-seenario  the  alloy  segregated  very  little.  Henee,  this 
set  of  DS-simulations  did  not  eonclusively  indieate  an  optimum  mesh. 

To  see  the  effeet  of  the  mesh  spaeing  on  simulated  maerosegregation  in  VAR- 
ingots,  a  refined  mesh  spaeing  was  employed.  The  dimension  of  the  fully-grown  domain 
is  43  cm  in  width  and  165  cm  in  height.  For  the  simulations  of  Figs.  4-6,  there  were  30 
elements  across  the  width  and  76  elements  along  the  height  of  the  eomplete  ingot,  so  each 
element  in  the  mesh  was  14.3  mm  x  21.7  mm.  In  this  seetion,  the  mesh  is  60  elements 
across  the  width  and  152  elements  along  the  height,  so  eaeh  element  in  the  mesh  is  7.2 
mm  X  10.8  mm. 
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1.6H 


Fig.  2:  Velocities  calculated  from  MULTIA:  (a)  after  6000  s  and  (b)  after  8700  s.  The 
extent  of  the  mushy  zone  is  also  shown.  Dimensions  are  in  m. 
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Fig.  3:  Macrosegregation  of  oxygen  calculated  from  MULTIA  after  8700  s  of  ingot 
solidification:  (a)  primary  arm  spacing  is  800  pm  and  (b)  primary  arm  spacing  is  300 
pm.  The  extent  of  the  mushy  zone  is  also  shown.  Dimensions  are  in  m. 
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Fig.  4:  Macrosegregation  of  oxygen  calculated  from  MULTIA  after  8700  s  of  ingot 
solidification  with  a  primary  arm  spacing  of  300  pm:  (a)  element  size  in  mesh  is  14.3  mm 
X  21.7  mm  and  (b)  7.2  mm  x  10.8  mm.  The  extent  of  the  mushy  zone  is  also  shown. 
Dimensions  are  in  m. 
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Fig.  5:  Macrosegregation  of  oxygen  calculated  from  MULTIA  after  16000  s  of  ingot 
solidification  with  a  primary  arm  spacing  of  300  pm:  (a)  element  size  in  mesh  is  14.3  mm 
X  21.7  mm  and  (b)  is  7.2  mm  x  10.8  mm.  Dimensions  are  in  m. 
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Fig.  6:  Velocities  calculated  from  MULTIA  with  a  primary  arm  spacing  of  300  pm  and 
element  size  of  7.2  mm  x  10.8  mm:  (a)  after  2000  s  and  (b)  after  4000  s.  The  extent  of 
the  mushy  zone  is  also  shown.  Dimensions  are  in  m. 
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Fig.  7:  Velocities  calculated  from  MULTIA  with  a  primary  arm  spacing  of  300  pm  and 
elements  of  7.2  mm  x  10.8  mm:  (a)  after  6000  s  and  (b)  after  8700  s.  The  extent  of  the 
mushy  zone  is  also  shown.  Dimensions  are  in  m. 
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Figures  4(a)  and  4(b)  show  the  macrosegregation  of  oxygen  after  8700  s  in  an 
ingot  with  a  primary  dendrite  arm  spacing  of  300  pm  used  for  predicting  the  permeability 
in  the  mushy  zone.  Figure  4(a)  is  the  result  simulated  with  the  larger  mesh  spacing  (14.3 
mm  X  21.7  mm),  and  Fig.  4(b)  is  with  the  smaller  mesh  spacing  (7.2  mm  x  10.8  mm). 
Qualitatively,  the  segregation  patterns  are  similar  except  more  details  of  the  segregation 
are  revealed  when  the  finer  elements  are  used.  Specifically  an  added  band  of  positive 
segregation  of  oxygen  exists  at  the  mid- width  of  the  ingot  for  the  simulation  with  the 
finer  mesh.  The  simulation  with  the  larger  elements  shows  somewhat  more  segregation, 
but  the  segregated  bands  are  narrower. 

Simulations  of  the  segregation  in  the  completely  solidified  ingots  are  given  in 
Figs.  5(a)  and  5(b).  Segregations  in  the  completely  solidified  ingots  are  similar  to  the 
results  shown  in  Fig.  4(a)  and  4(b),  respectively,  which  indicates  that  the  pattern  of 
segregation  was  essentially  developed  at  the  earlier  time.  Near  the  right  side  of  each 
ingot,  there  is  a  band  depleted  in  oxygen  and  enriched  in  iron.  In  the  previous  section,  it 
is  reported  that  the  vertical  bands  are  possibly  weak  freckles,  but  the  DS-set  of 
simulation-experiments  offers  evidence  that  the  Ti-6A1-4V  alloy  has  no  tendency  for 
forming  vertical  freckles.  Figures  6(a)  and  (b)  and  Figures  7(a)  and  7(b)  show  the 
calculated  velocities  in  the  melt  pool;  these  figures  represent  output  after  2000  s,  4000  s, 
6000  s  and  8700  s  of  solidification  of  the  ingot.  Because  of  the  Lorentz  forces,  the  melt- 
pool  sweeps  upward  along  the  leading  part  of  the  mushy  zone  and  carries  oxygen- 
depleted  melt  toward  the  upper  right  comer,  where  there  are  weak  secondary  convection 
cells.  Hence,  a  depletion  of  oxygen  in  the  right  side  of  the  ingot  results,  but  the  exact 
mechanism  making  the  segregation-band  near  the  right  side  of  ingot  is  not  known. 
Perhaps  the  mechanism  is  related  to  the  fact  that  the  melt  pool  is  almost  stagnant  where 
the  vertical  oxygen-depleted  band  intersects  the  leading  part  of  the  mushy  zone. 

IV.  Status  of  VAR-Simulations 

In  this  report  it  is  demonstrated  that  MULTIA  can  be  used  to  predict  the 
formation  of  macrosegregation  in  ingots.  In  its  present  form,  however,  MULTIA  is  an 
"in-house"  research  code  without  requisite  “front-end”  and  “outpuf  ’  set-ups  for 
engineering  applications.  There  are  computer  codes  for  simulating  VAR-ingots,  but  they 
are  not  used  by  USA  VAR-ingot  producers  for  predicting  and  controlling 
macrosegregation.  The  codes  are  called  BAR,  Calcosoft,  SOLID,  and  COMPACT. 

BAR  is  a  proprietary  code  of  the  Specialty  Metals  Processing  Consortium 
(SMPC).  The  industrial  members  are  companies  who  produce  VAR-  and  ESR-ingots  of 
titanium  alloys,  superalloys,  and  specialty  steels.  BAR  solves  the  energy  equation,  the 
momentum  equation  with  a  simple  turbulence  model,  and  the  Lorentz  forces,  but  does  not 
calculate  the  redistribution  and  transport  of  the  alloy-  and  impurity-elements.  It  does  not 
calculate  macrosegration,  and  it  is  used  exclusively  for  VAR  processing.  Simulation  runs 
using  BAR  are  very  fast,  so  it  is  a  candidate  for  real-time  control  systems. 

SOLID  is  software  developed  at  the  Laboratory  of  Science  and  Engineering  of 
Materials  and  Metallurgy  in  the  School  of  Mines  in  Nancy,  France.  This  code  can  be  used 
exclusively  by  companies  that  are  contributing  members  to  the  Solidification  Group  at 
the  School  of  Mines.  SOLID  simulates  the  solidification  of  ingots,  including  VAR-ingots, 
and  does  calculate  macrosegregation.  Among  the  SMPC  members.  Timet  is  a  member 
company  of  both  the  School  of  Mines  in  Nancy  and  the  SMPC  and  uses  SOLID. 
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According  to  a  user:  Calcosoft  is  licensed  by  the  European-based  software 
company  known  as  the  ESI  group.  "CALCOSOFT-2D  is  a  promising  tool  for  the 
modeling  and  simulation  of  remelting  processes  like  Vacuum  Arc  Remelting  (VAR). 
Based  on  our  current  simulation  results  and  on  dedicated  software  developments  in  the 
frame  of  a  joint  project  with  Calcom  ESI,  we  will  be  able  to  improve  the  quality  of  our 
remelted  ingots  by  an  optimization  of  the  melt  rate,  the  fdl  ratio,  the  cooling  system  and 
the  current  and  voltage  definitions."  [Christof  Sommitsch,  Project  Leader,  Bohler 
Edelstahl  GmbH  &  Co  KG,  Austria  .]  According  to  the  Calcosoft-ESI  website,  a  module 
for  calculating  macrosegregation  is  available. 

COMPACT  is  a  product  of  Innovative  Research,  Inc.  (IRI).  Under  the 
COMPACT  umbrella  customers  choose  a  Customized  Package.  SMPC  included  an 
evaluation  of  COMPACT  for  ESR-ingots  as  part  of  a  year  long  program.  The  Customized 
Package  delivered  by  IRI  was  suitable  for  “steady-state”  ingot  solidification.  As  of  2005 
a  version  for  analyzing  ingots  that  solidify  under  transient  conditions  had  not  been 
released  by  IRI.  More  importantly  COMPACT  does  not  calculate  macrosegregation.  No 
USA  company  is  using  COMPACT  for  VAR-processing.  A  version  of  the  transient  VAR 
model  is  under  development  so  that  situation  may  change.  An  SMPC  member  believes 
there  are  three  USA  licensees  for  ESR  processing. 

Under  the  aegis  of  SMPC,  Purdue  investigators  have  used  PHYSIC  A  of  Physica 
Ltd.,  U.K.  to  simulate  macrosegregtion  in  binary  alloys  of  Ti  under  VAR  conditions. 
PHYSICA  is  not  typically  used  for  ingot  solidification,  but  through  an  arrangement  that 
allowed  the  Purdue  investigators  access  to  the  object  code,  they  were  able  to  design 
solidification-modules  that  were  combined  with  PHYSICA  to  effect  the  simulations.  The 
Purdue/Physica  arrangement  was  suitable  for  carrying  out  research,  but  neither  SMPC, 
Purdue  nor  Physica  Ltd.  has  commercialized  specialized  software  for  VAR-ingot 
producers.  None  of  the  SMPC  members  uses  PHYSICA  for  simulating  VAR-  or  ESR- 
ingots. 

As  of  this  writing,  it  appears  that  only  CALCOSOFT-2D  could  be  used  by  the 
VAR-ingot  producers  in  the  USA  for  calculating  macrosegregation.  Timet  of  Henderson, 
Nevada  is  an  exception  because  Timet  is  a  member  company  of  the  School  of  Mines  in 
Nancy,  France.  The  members  of  SMPC  use  BAR,  but  of  course  cannot  do 
macrosegregation  calculations  with  it.  MULTIA  was  originally  designed  to  simulate  the 
formation  of  macrosegregation  (including  freckles)  in  multi-component  alloys  that  are 
directionally  solidified.  As  mentioned  above,  in  order  for  MULTIA  to  be  used  as  an 
engineering  tool  for  simulating  Ti-alloy  ingots  there  are  several  to  many  features  that 
should  be  added.  Appendix  B  summarizes  enhancements  to  make  MULTIA  into  an 
engineering  tool. 
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Appendix  A  -Continuum  Model  of  Dendritic  Solidification 

(Appendix  A  was  prepared  by  D.R.  Poirier,  Department  of  Materials  Scienee  and 
Engineering,  The  University  of  Arizona  and  S.D.  Felieelli,  Mechanieal  Engineering 
Department,  Mississippi  State  University.) 

A.l  Introduction 

This  progress  report  details  the  eontinuum  model  that  is  the  basis  of  the 
solidifieation  simulator  at  The  University  of  Arizona,  which  is  known  as  MULTIA.  The 
underlying  model  comprises  the  set  of  transport  equations:  continuity,  momentum, 
energy  and  solute  conservation.  The  mushy  zone  is  treated  as  a  porous  medium  in  the 
model.  The  numerical  method  that  is  used  to  implement  the  model  into  MULTIA  is 
summarized.  MULTIA  is  a  finite-element  code  that  is  used  to  simulate 
macrosegregation,  especially  the  formation  of  freckles,  and  porosity- formation  in 
multicomponent  cast-alloys.  In  this  report,  the  derivations  of  the  conservation  equations 
and  the  main  features  of  the  numerical  methods  and  algorithms  that  underlie  MULTIA 
are  summarized. 

A.2  Continuity  Equation 

The  following  assumptions  are  invoked: 

(i)  the  solid  is  stationary  and  only  the  liquid  convects; 

(ii)  the  densities  of  the  liquid  and  solid  differ  and  each  is  constant. 

Assumption  (i)  leads  to 

d 

^Xpsgs+ PLgJ  +  ^  ■  Pl^  =  0  (Al) 

ot 
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where  t  is  time  (s);  and  are  the  densities  of  the  solid  and  liquid  (kg  m'^) 
respectively;  and  the  volume  fraction  of  solid  and  liquid,  respectively;  and  u 

is  the  superficial  velocity  (m  s  ’).  Equation  (1)  can  be  found  in  Ganesan  and  Poirier  [Al], 
In  a  porous  medium  {i.e.,  the  mushy  zone)  the  superficial  velocity  is  the  volume  flow  rate 
in  a  direction  divided  by  the  unit  area  perpendicular  to  the  flow.  The  unit  area  contains 
both  the  solid  and  liquid.  In  the  liquid  pool,  «  is  simply  the  velocity,  and  in  the 
completely  solidified  part  of  the  ingot  «  =  0. 

Since  the  densities  of  the  solid  and  liquid  are  each  constant,  Eq.  (Al)  simplifies  to 

Ps  -^  +  Pl  +  Pl  V  • «  =  0 
®  dt  dt 

but  dg^/dt  =  -dg^/dt ,  so  that 

,  _V.„ 

Pl 

We  define  the  solidification  shrinkage  as 


A  = 


Ps  -  Pl 
Pl. 


so  the  final  form  of  the  continuity  equation  becomes 


V-«  =  /3^ 
dt 


where  (p  =  gi^. 


(A2) 
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A.3  Momentum  Equation 

Ganesan  and  Poirier  [Al]  used  the  teehnique  called  volume  averaging  to  derive 
the  momentum  equation  that  applies  to  the  mushy  zone,  which  is  treated  as  a  porous 
medium.  The  following  assumptions  are  invoked: 

(iii)  the  exchange  of  momentum  because  of  solidification  shrinkage  is  neglected; 

(iv)  the  average  momentum-dispersion  is  neglected; 

(v)  the  Newtonian  constitutive  relation  between  stress  and  strain  rates  is  assumed;  and  it 
is  assumed  that  first  order  resistance  to  flow  is  adequate  to  describe  flow  (this 
amounts  to  Darcy  flow). 

Under  these  assumptions  and  according  to  Ganesan  and  Poirier  [Al],  the  momentum 
equation  is 

^  =  -S^yp  +  +  V-7V  -//g/[K]-‘v  (A3) 

_dt 

where  v  is  the  intrinsic  velocity  of  the  interdendritic  liquid,  P  is  the  pressure  (Pa),  g  is 
gravitational  acceleration  (m  s'^),  //  is  the  viscosity  of  the  interdendritic  liquid  (Pa  s),  K 
is  the  permeability  tensor  (m^),  and  N  arises  from  volume  averaging  of  the  deviatoric 
part  of  the  stress  tensor  of  the  interdendritic  liquid. 

The  intrinsic  velocity  is  the  volume-averaged  velocity  of  the  interdendritic  liquid; 
it  is  related  to  the  superficial  velocity  by 

u 

V  =  — 

<!> 

The  fourth  term  on  the  right  side  of  Eq.  (A3)  accounts  for  the  resistance  to  flow  offered 
by  the  dendritic  network  within  the  mushy  zone.  The  permeability  K  is  a  conductivity 
for  flow  in  the  mushy  zone,  which  depends  on  the  volume  fraction  of  liquid  and  a  length 
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scale  of  the  dendritic  solid.  Where  (/)=0  (/.e.,  all  solid),  K  =  0,  v  =  «  =  0  and  the 
momentum  equation  does  not  apply.  Where  ^=1  {i.e.,  the  liquid  pool),  all  components 
of  K  are  oo,  and  Eq.  (3)  reduces  to  the  usual  Navier-Stokes  equation  for  a  liquid  [A2]. 
Thus  Eq.  (A3)  allows  for  a  natural  transition  in  momentum  transport  from  the  mushy- 
zone  to  the  liquid  pool.  Since  (j)  is  treated  as  an  independent  variable,  the  position  of  the 
leading  front  of  the  mushy  zone  is  calculated  as  part  of  the  solution.  Furthermore,  the 
liquid  at  the  leading  front  of  the  mushy  zone  is  not  assumed  to  be  of  the  nominal 
composition  and  at  the  liquidus  temperature  of  the  nominal  alloy. 

To  include  the  effect  of  shrinkage  we  start  with  the  stress  tensor  in  Eq.  (A3).  It 
can  be  written  as 


V-iV 


//V- 


v«  +  {v«}^ 


2 

3 


{V  • 


(A4) 


where  the  second  term  on  the  right  side  is  the  transpose  of  the  matrix  Vh  ,  and  is  the 
second  order  unit  tensor. 

For  convenience,  each  term  in  Eq.  (A4)  is  expanded  in  two  dimensions  {x,y)  with 
velocity  components  {u,  v): 


V« 


{Vuf 


du 

dv 

dx 

dx 

du 

dv 

dy 

du 

du 

dx 

dy 

dv 

dv 

dx 

and 
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(2)  _ 


du  dv 
dx  dy 


du  dv 
dx  dy 


After  summing  the  three  matrices,  we  get 


V-N  =  // 


^  ]  d  ^ 

(V-«) 

3  dx 


,  1  d 

V^v  +  (V-«) 

3  dy 


or  generalizing  to  three-dimensions: 

1  d 


V-N  =  jLi 

+  M 

+  M 


+ 


3  dx 

1  a 


{V-u) 


V^v  +  (V-«) 

3  dy 

'  .  1  a 

+ - (V-«) 

3  az 


(x  -  component) 
{y  -  component) 
(z  -  component) 


(A5) 


where  the  components  of  u  are  (w,  v,  w).  Notice  that  each  term  contains  the  divergence 
of  «;  continuity,  Eq.  (A2),  is  substituted  so  that 


V-N  =  ju 


3  dx  dt 


+  A 


V^v  + 

3  dy  dt 


+  A 


3  az  dt 


or 


V-N 


//  V^u  + 


J 

3  V  y 


(A6) 


Equation  (A6)  is  substituted  into  Eq.  (A3)  to  get  an  equation  for  the  conservation  of 
momentum  with  shrinkage;  the  notation  change  of  (p  for  is  also  made  and  u  =  (/)v  is 


substituted.  The  result  is 
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du  I 

—  +  «-V 
dt 


r . 


-^VP  +  (/) 
Pl 


^  P  1  V72  Tli 

—  s  +  vV  «  +  — 

PlJ  3 


v^[K]  ‘« 
(A7) 


Equation  (A7)  is  the  same  as  Eq.  (2)  in  Felicelli  et  al.  [A3]. 

The  only  remaining  detail  in  the  momentum  equation  is  the  Boussinesq 
approximation.  It  is  assumed  that  the  density  of  the  liquid  is  constant  at  its  reference 
value,  yC>L ,  in  all  terms  except  in  the  body  force  term  {i.e.,  the  gravity  term).  For  that  we 


use 


where 


and 


_P_ 

Pi. 


1  +  PrlT-T,)  +  'Z/S’iC'-Ci) 

j  =  l 


Pt 


1  dp 


1  dp 

■ 


(A8) 


In  Eq.  (A8)  N  is  the  number  of  components,  j  refers  to  a  component,  and  Tj^  and  C{ 
refer  to  the  reference  temperature  and  concentration,  respectively.  The  reference  state 
has  the  coordinates: 

T^,ClCl,Cl . 


where  is  the  liquidus  of  the  alloy  with  its  nominal  composition  and  the  Q  ’s  refer  to 
the  concentrations  of  the  elements  1,  2,  3 . in  the  nominal  alloy-chemistry. 
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A.4  Solute  Conservation  Equation 

Assumption  (i)  still  applies.  Other  assumptions  at  the  local  level  are  explained 
below.  Starting  with  Eq.  (15)  from  Poirier  et  al.  [A4],  each  component  in  the  alloy  or 
dissolved  gas-forming  element  follows 

V.y  +  V.(p,Q«)  +  |-(^)  =  0  (A9) 

(At  this  point,  the  superscript  j  has  been  omitted  to  simplify  notation.  When  it  is 
necessary  to  consider  all  of  the  components  together,  the  superscript  is  included.)  In  Eq. 

(A9),  j  is  the  diffusion  flux  in  kg  (solute)  m'^  s  ’,  pC  is  the  concentration  in  kg  (solute) 
m'  ,  and  the  second  term  is  the  advection  of  solute  because  the  liquid  has  velocity  u. 
Equation  (A9)  applies  to  the  mixture  (z.e.,  mushy  zone)  and  the  all-liquid  zone.  With  u  = 
0  it  also  applies  to  the  all-solid  zone. 

For  our  assumptions  of  constant  and  constant  (3 ,  the  advection  term  becomes 


(AlO) 


(All) 


(A12) 
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where  C  =  pCj and  Cg  is  the  average  concentration  of  the  solute  in  the  solid 

(kg  (solute)  kg'^  or  wt.%).  The  model  considers  two  types  of  elements:  those  in  which  it 
is  assumed  there  is  no  diffusion  in  the  solid  and  those  in  which  there  is  complete 
diffusion  in  the  solid. 

For  the  alloy  elements  with  no  diffusion  in  the  solid,  the  solid  that  forms  from  the 
liquid  is  k  with  k  defined  as  the  equilibrium  partition  ratio.  Then  the  average 
concentration  of  an  element  in  the  solid  (wt.%)  is 

Cs  =  <A13) 

When  the  model  is  incorporated  into  MULTIA,  it  is  necessary  to  evaluate  and  store  the 
integral  at  every  node  in  the  finite-element  mesh. 

For  interstitial  elements  with  complete  diffusion  in  the  solid: 

Cg  =  Q  =  kC^  (A14) 

Returning  to  Eq.  (All),  we  see  that  the  flux  for  diffusion  in  the  mixture  is  needed.  This 
is  written  as 

j  =  -p^D^VC^  -  p,D,VC,  (A15) 


where  and  Z)g  are  effective  diffusivities  in  the  liquid  and  solid.  These  are 


= 


A  = 


24 


where  and  Dg  are  the  diffusion  coefficients  and  and  rg  are  the  tortuosities  to 
account  for  the  morphology  and  connectivity  of  the  phases.  Lacking  data  on  the 
tortuosities,  we  assume  that  =  Tg  =  1  (Assumption  vi). 

For  most  alloy  elements  {i.e.,  substitutional  elements),  »  Dg  so  Eq.  (A  15) 
reduces  to 

j  =  -p^(l>D^WC^  (A16) 

On  the  other  hand,  some  elements  diffuse  rapidly  in  the  solid  {i.e.,  interstitial  elements 
such  as  N,  H,  O  and  C  in  nickel  alloys,  titanium  alloys,  and  ferrous  alloys  and  H  in 
aluminum  alloys).  These  elements  maintain  equilibrium  during  partitioning;  thus 

VCg  =  kVQ 

so  that  Eq.  (A15)  becomes 

j  =  +  kps{l-^)D^]VC^ 

or 

J  =  (A17) 

Thus,  the  effective  diffusivity  of  the  mixture  is 

^mix  ~  ~  [  Pl  ^  ^  Ps  (1  ~  ] 

p 

and  p  is  the  mixture  density: 

P  =  <t>  Pi.  ^ 

Equation  (A17)  is  the  same  as  Eq.  (1 1)  in  Ref  [A3]. 
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A.5  Conservation  of  Energy 

Equation  (1)  in  Poirier  et  al.  [A4]  can  be  written  as 


^(pH)  =  v.(/rvr)  -  y-{p^H^u) 

ot 

where 


(A18) 


pH  =  p,{\-<j>)H,  +  p^c^H^  (A19) 

K  is  the  mixture  thermal  conductivity,  u  is  the  superficial  velocity  of  the  liquid,  and 
and  are  the  intensive  enthalpies  of  the  solid  and  liquid  (J  kg  ').  The  enthalpies  are 

expressed  in  terms  of  the  respective  enthalpies  at  a  reference  temperature,  ,  where  the 
heat  of  fusion  is  given.  Hence 

H,  =  +  c,(T-T^)  (A20) 

//l  =  +  c^{T-T^)  (A21) 

and 

L  =  Hi  +  Hi  {All) 

where  Cg  and  are  the  specific  heat  capacities  (J  kg''  K'')  and  L  is  the  latent  heat  of 
fusion  (J  kg’'). 

The  time  derivative  of  the  volumetric  enthalpy  in  Eq.  (A18)  is  expanded  as 
follows: 


dH^ 


—  [ph]  =  -PsH^  +  Pl<^ 

dr  ’  ^  ^  dt  ^  ^  dt  ®  dt  ^  dt 


(A23) 


The  advective  term  in  Eq.  (A18)  is  expanded  while  taking  p-^  =  constant: 
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+  up^VH^ 


=  +  Pl^QVJ  (A24) 

Now  we  combine  Eqs.  (A18),  (A23)  and  (A24): 

[p^{\-^)c,+p^<l>c^]^  =  V -{kVT)- p^c^u-VT  -  p,[l  +  {c^-c,){T -T^)]^ 

(A25) 

Numerically,  it  is  important  to  express  the  term  involving  d(f>ldt  as  a  function  of 
dTidt . 

The  latent  heat  term  in  Eq.  (A25)  must  be  treated  implicity  in  order  for  the 
numerical  algorithm  in  MULTIA  to  be  stable  [A5-A7].  To  make  this  term  implicit  for  a 
multicomponent  alloy  requires  that  the  temperature  of  the  mixture  is  related  to  the 
composition  of  the  interdendritic  liquid  [A8,A9]. 

Consider  first  a  solute  that  diffuses  completely  in  the  local  solid,  so  that 
Q  =  A:  Cl  in  Eq.  (A12).  Equation  (A12)  can  then  be  written 

C'  =  [k^  + 

where  ={\  +  /3)  and  j  refers  to  a  particular  element.  Hence, 

1  dC'  _  ^ 

K^+(l-K^')^  dt  ~  dt  Rj +(l-K^')^  dt 


(A26) 


(A27) 


The  temperature  in  the  mushy  zone  is  the  liquidus  temperature  of  the  interdendritic 
liquid,  so  that 


dt 


Z 


m 


dCj 

dt 


(A28) 
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where  =dTj dCl  represents  the  change  in  the  temperature  as  €[  changes  and  N 

refers  to  the  number  of  elements  in  the  alloy. 

N 

For  a  reason  that  will  become  apparent,  Eq.  (A27)  is  multiplied  by  . 

y  =  i 


N 


dC 


fr',  +{\-K^)(l)  dt 


=  y  +  ]y  d(/> 


(A29) 


By  Eq.  (A28),  the  first  term  on  the  right  side  is  simply  dTjdt .  To  simplify  the  notation, 


we  define: 


M-'  = 


m 


and 


Finally,  Eq.  (A27)  yields: 


=  (l-K^)M^Cl. 


iM  =  If  M 

dt  "  ^ 


A 

A  7=1 


dc 

dt 


J_^ 

dt 


(A30) 


At  this  point  we  eliminate  d^jdt  in  Eq.  (A25)  by  substituting  Eq.  (A30);  before 
doing  so,  we  divide  Eq.  (A25)  by  .  These  two  steps  result  in 


,dT 


[r  +  (l  +  j0)(l-^)  -  aV^T  =  -u  VT  + 


(A31) 


where  /  =  c^  I c^',  is  assumed  to  be  constant;  a  =  Kj and  A^  is  defined  as 


A,  =  (\  +  P)[{Llc^)  +  {\-y){T-T^)-\. 


By  combining  Eqs.  (A30)  and  (A31),  we  obtain  the  energy  equation  without  , 
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(A32) 


[y{\  +  P)(\-cI>)  +  cI>-A]—  -  =  u-VT  -  ^ 

ut  7=1 

where  A  =  A^j A^  in  the  mushy  zone  and  ^  =  0  in  the  all-liquid  and  all- solid  zones. 
Recall,  however,  that  the  parameter  (and  hence  A^ )  depends  on  the  assumption  that 
the  solute  elements  diffuse  completely  in  the  local  solid. 

It  remains  to  determine  the  energy  equation  when  there  is  no  diffusion  in  the  local 
solid.  For  this  case  we  return  to  the  equation  just  above  Eq.  (A12)  and  expand  it  as 
follows: 

|■(^)  =  (A33) 

Equation  (A33)  replaces  Eq.  (A27)  in  the  derivative  of  the  final  form  of  the  energy 
equation.  The  result  is  that  Eq.  (A32)  also  applies  for  no  diffusion  in  the  solid,  but 
(and  hence  A^ )  takes  on  different  values: 

=  —  (A34) 

Since  is  included  in  the  definition  of  A^ ,  then  it  too  is  affected.  Equation  (A32)  was 
first  presented  in  Ref  [A3]  but  without  its  derivation. 

A.6  Numerical  Method 

The  governing  equations  are  solved  with  a  finite  element  discretization,  using 
bilinear  isoparametric  elements,  and  a  finite  difference  scheme  to  discretize  in  time.  The 
algorithm  uses  a  Petrov-Galerkin  formulation  for  stabilization  of  the  convective  terms, 
and  the  penalty  method  to  impose  incompressibility.  The  algorithm  solves  the 
conservation  equations  written  in  terms  of  nondimensional  variables;  many  details  are  in 
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Refs.  [A7]  and  [A9]  and  are  not  repeated  here.  However,  differences  in  the  numerical 
expressions  and  calculations  procedures,  arising  because  of  the  inclusion  of  interstitial 
elements  and  the  effect  of  shrinkage,  need  to  be  mentioned.  Also,  in  addition  to  the 
formal  procedures  involved  in  the  finite  element  method,  a  successful  numerical 
implementation  of  the  solidification  model  requires  special  consideration  and 
organization  of  the  calculations. 


A.6.1  Computation  of  the  Solute  Concentrations  in  the  Liquid  and  Solid 

The  average  solute  concentration  C  ^  is  computed  from  the  conservation 
equation,  Eq.  (A1 1).  Using  this  value,  we  calculate  the  concentrations  in  the  phases 
following  the  procedure  described  in  Ref  [A9].  Some  modifications  need  to  be  done, 
however,  to  allow  for  shrinkage  and  solutes  that  have  high  diffusion  in  the  solid  (like  gas¬ 
forming  elements).  Defining 

P  =  \\k^Cid<j>  =  {\-<j>)Ci  (A35) 


and  proceeding  as  in  Ref  [A9],  we  arrive  at  the  following  expressions  for  calculating  the 
phase  concentrations  of  solutes  that  have  negligible  diffusion  in  the  solid: 

(1)  Solidification 


rJ,  «+l  _ 


C’"  +  (Cl’"  +  Cl’"^f(f 


c 


j,  n+l 


c 


j,  n+\ 


(\  +  +\k^ Cl’"  (f 


(A36) 

(A37) 


(2)  Remelting 


(A3  8) 
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n+\ 


(A39) 


<!> 


n+\ 


where  the  subscript  n  indicates  the  time  level  f ,  and  is  a  function  that  stores  the 
histories  of  solidification  for  each  alloy  component. 

For  the  case  of  solutes  that  have  high  diffusivities  in  the  solid,  the  following 
expressions  are  obtained  for  both  solidification  and  remelting  situations: 


(A40) 


Q  7,  n+\ 


(A41) 


where  =i\  +  fi)  .  For  simplicity  of  notation,  the  partition  ratios  are  assumed 
constant  in  Eqs.  (A36)  through  (A41);  however,  the  model  allows  them  to  be  functions  of 
concentration,  and  hence  of  time. 


A.6.2  Computation  of  the  Volume  Fraction  of Liquid 

Here  we  depart  from  the  method  used  in  Ref  [A9],  which  expressed  the  liquidus 
temperature  in  a  pseudo-lineal  form  and  hence  required  that,  for  general  nonlinear 
liquidus  functions,  modified  slope  coefficients  m^  be  calculated,  in  addition  to  the 
regular  coefficients  .  This  is  not  needed  in  the  current  approach,  and  arbitrary  liquidus 
functions  can  be  easily  accommodated  in  the  model.  Since  the  concentrations  of  the  gas¬ 
forming  elements  are  very  small,  we  neglect  their  effects  on  the  liquidus  temperature  of 
the  alloy.  The  concentration  of  solute  j  in  the  liquid  is 

^  C  ^  -  (1  +  A) 

'  <!> 

and  substituting  this  into  the  liquidus  function,  we  have 
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Tu,  =  fiCi)  =  f 


^  -  (1  +  /?)  P  ^ 

</> 


-  F{^) 


The  temperature  is  the  liquidus  temperature  that  is  a  function  of  the  composition  of  the 
interdendritic  liquid.  Hence, 

T  -  F{(l))  =  0  (A42) 

and  Eq.  (A42)  is  solved  using  the  secant  method.  At  time  level  and  for  every  node  in 
the  mushy  zone,  we  have  the  following  iteration  scheme: 

g'  =  T  -  Fif) 


=  S'  -  g 


For  iteration  /  =  0  : 


/  =  f-  g-‘  =  0 


-1 


where  S"  is  the  volume  fraction  of  liquid  at  the  previous  time  level  f 


A.6.3  Calculation  of  Pressure 

To  obtain  the  velocity  field,  the  weighted  residual  forms  of  Eqs.  (A2)  and  (A7) 
(in  nondimensional  form)  are  discretized  using  bilinear  isoparametric  elements  and  a 
penalty  method  to  impose  incompressibility.  The  use  of  the  penalty  method,  however, 
eliminates  the  pressure  as  an  unknown,  which  needs  to  be  recovered  through  a  post¬ 
processing  step.  In  addition,  because  the  continuity  equation  ( Eq.  (A2) )  is  not 
homogeneous,  the  implementation  of  the  penalty  method  differs  from  the  standard 
approach  used  in  the  incompressible  Navier-Stokes  equations.  To  attack  this  problem, 
we  use  the  method  proposed  by  McBride  et  al.  [A  10].  They  showed  that  discretizing  the 


32 


continuity  and  momentum  equations  with  a  mixed  formulation,  consisting  of  a  bilinear 
velocity  and  constant  element  pressure  at  the  element  level,  is  equivalent  to  the  use  of  a 
penalty  method  with  the  pressure  defined  by  the  relation  (in  nondimensional  form): 


P  =  Ps  -  — 

Re  I  at 


(A43) 


where  p  is  the  total  pressure,  is  the  static  pressure,  Rq  =  VH is  a  Reynolds 

number,  and  X  is  the  penalty  parameter.  The  penalty  parameter  must  be  a  large  number 
(e.g.,  X  =  10’°)  in  order  to  impose  incompressibility  effectively.  This  number  may  need 
to  be  even  larger  when  the  range  of  permeability  in  the  mushy  zone  varies  greatly,  and 
the  Darcy  term  becomes  very  large  for  small  fractions  of  liquid  in  the  mushy  zone. 


A.6.4  Program  Algorithm 

The  governing  equations  are  solved  sequentially,  and  a  partial  iteration  is 
performed  within  each  time-step  to  obtain  convergence.  The  iteration  does  not  include 
the  momentum  equations,  which  are  solved  only  once  at  the  beginning  of  the  time-step. 
At  time  f  all  conditions  are  known.  To  advance  to  time  =  t"  -h  At  and  whenever  an 
equation  is  solved  for  one  of  the  dependent  variables,  the  latest  available  values  for  all 
other  variables  are  used.  The  following  scheme  is  used  in  the  calculations: 

(1)  Solve  for  velocity  and  pressure,  from  the  nondimensional  momentum  equation  and 
Eq.  (A43). 

(2)  Solve  for  temperature  from  the  nondimensional  energy  equation. 

(3)  Compute  fraction  of  liquid  from  Eq.  (A42). 

(4)  Compute  phase  concentrations  from  Eqs.  (A36)  through  (A41). 

(5)  Solve  for  mixture  concentrations  from  the  nondimensional  form  of  Eq.  (A1 1). 
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Step  1  is  solved  onee  per  time  level,  while  steps  2  through  5  are  iterated  within  a 
time  level  until  convergence.  Further  details  on  the  algorithm  can  be  found  in  Refs.  [A7] 
and  [A9],  including  treatment  of  remelting,  solidification  at  the  eutectic  temperature,  and 
selection  of  time-steps  and  meshing  strategies. 
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Appendix  B  -  Enhancements  to  MULTIA  for  Simulating  Macrosegregation  in  VAR- 
Ingots 

1 .  Building  an  Ingot.  In  the  course  of  this  program,  the  capability  of  adding  material  to 
the  calculation  domains  the  ingot  is  formed  from  the  consumable  electrode  was 
implemented 

2.  Cylindrical  Coordinates.  In  this  program,  the  MULTIA  demonstrations  were  for  an 
ingot  shaped  as  a  slab.  Rectangular  coordinates  were  used.  For  VAR  applications, 
however,  cylindrical  coordinates  are  applicable.  MULTIA  can  be  run  with  cylindrical 
coordinates,  but  the  matrices  are  ill-conditioned  and  the  execution  time  is  excessive. 

There  are  better  techniques  for  solving  the  matrices  that  result  when  cylindrical 
coordinates  are  used,  and  these  should  be  pursued. 

3.  Motion  of  Equiaxed  Grains.  Presently  MULTIA  assumes  that  the  solid  in  the  mushy 
zone  is  stationary.  The  solidifying  equiaxed  grains  in  VAR-ingots,  however,  are  free  to 
convect  with  the  liquid  until  the  fraction  of  solid  achieves  a  critical  concentration,  above 
which  the  grains  develop  a  stationary  network  of  the  dendritic  grains.  The  motion  of  the 
solid  can  contribute  to  the  overall  segregation,  so  incorporating  this  into  MULTIA  would 
improve  macrosegregation  predictions. 

4.  Lorentz  Forces.  In  this  program,  the  Lorentz  forces  were  transferred  from  BAR  and 
added  as  body  forces  in  MULTIA.  It  would  be  relatively  easy  to  incorporate  the  electro¬ 
magnetics  to  make  MULTIA  independent  in  this  regard,  once  correct  boundary 
conditions  for  current  distribution  are  established. 

5.  Thermal  Boundary  Conditions.  In  this  program,  the  boundary  temperatures  were 
transferred  from  BAR  to  MULTIA.  Boundary  conditions  in  the  VAR-process  are 
complex  and  should  be  thoroughly  assessed.  The  major  complexities  include  heat  transfer 
coefficients  that  vary  with  time  and  location  as  a  gap  develops,  inclusion  of  the  bottom 
stool  (starting  block)  in  the  calculation-domain,  energy-  and  material-transfer  from  the 
consumable  electrode,  proper  boundary  conditions  to  capture  radiation  exchange  among 
the  melt-surface,  electrode  and  crucible  wall,  etc. 

6.  Conformable  Meshing.  A  major  restriction  is  that  the  mesh  should  be  refined  enough 
so  that  between  the  liquidus  and  the  solidus  of  the  mushy  zone  there  are  several  elements 
(say  5  to  10)  at  a  given  elevation  in  the  ingot.  The  mushy  zones  in  Ti-ingots  are  rather 
thin  which  means  that  in  a  uniform  mesh  the  total  number  of  elements  in  the 
computational  domain  can  be  very  large.  Hence,  conformable  meshing  to  accommodate 
small  elements  in  the  mushy  zone  and  larger  elements  in  the  liquid  pool  and  solidified 
portion  of  the  ingot  should  be  adapted  to  MULTIA  to  keep  computational  times  within 
reason. 

7.  Turbulence.  MULTIA  assumes  laminar  flow  in  the  liquid  pool.  Turbulence  models 
should  be  assessed:  the  one  presently  in  BAR  with  "eddy"  transport  properties;  simple 
estimates  of  mixing  lengths  for  Reynolds'  stresses;  the  often  applied  "A:-e"  model,  which 
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is  used  in  many  commercial  codes;  and  a  "large  eddy"  model.  These  are  listed  in  order  of 
sophistication  and  computational  time!).  The  last  in  the  list  has  been  used  to  simulate 
melt  flows  (  forced  convection )  in  continuous  casting  of  steel.  In  VAR-pools,  however, 
flows  are  driven  by  body  forces  (sum  of  the  buoyancy-  and  Lorentz-forces). 

8.  Parallelization.  Even  in  its  present  form,  MULTIA  runs  for  tens  of  hours  to  make  a 
simulation.  With  the  enhancements  outlined  herein,  it  is  expected  that  computational 
times  would  increase.  Hence,  MULTIA  should  be  set  up  for  multi-processor  systems. 
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Appendix  C  -  Formation  of  Microporosity  in  Thin  Wall  Investment  Castings  of 
IN718 

A  series  of  Progress  Reports  on  porosity  in  experimental  plate-castings  of  IN718, 
along  with  details  of  a  solidification  model,  were  submitted.  The  contract  support  from 
AFML  was  supplemented  with  an  unrestricted  grant  from  Honeywell  Engines  in  Phoenix, 
AZ.  The  progress  reports  are  listed  in  Table  Cl,  and  the  major  findings  are  summarized. 


Table  Cl.  Progress  Reports  under  Award  F336 15-02- 1-52 19,  “Formation  of  Micro¬ 
porosity  in  Thin  Wall  Investment  Castings  of  IN718.” 


Date  of  Report 

UA  Authors 

Subiectfs) 

09/13/2003 

Poirier  and  Sung 

Solidification  model 
and  simulations 

11/02/2003 

Lyman  and  Poirier 

Fluorescent  penetrant 
(FPI)  and  radiography 

01/09/2004 

Lyman  and  Poirier 

Porosity  assessment  via 
quantitative  metallography 

02/25/2004 

Lyman  and  Poirier 

Comparison  of  FPI  and 
quantitative  metallography 

04/12/2004 

Lyman,  Sung 
and  Poirier 

Distribution  of  porosity 

06/25/2004 

Lyman,  Sung 
and  Poirier 

Grain  size,  arm  spacings  and  porosity 

Experimental  castings  of  IN718  were  made  by  investment  casting  under  vacuum 
at  two  Howmet  Casting  facilities:  the  Howmet  Research  Center  in  Whitehall,  Michigan 
and  a  foundry  in  LaPorte,  Indiana.  In  Whitehall  experimental  plates  of  thickness  0.04, 
0.08  and  0.12  inches  were  cast  horizontally.  The  set  of  castings  made  in  LaPorte  was 
designed  by  Honeywell  Engines  as  part  of  an  in-house  program  on  assessing  the  efficacy 
of  thin  wall  castings.  In  LaPorte,  the  experimental  plates  of  thickness  0.04  and  0.055 
inches  were  cast  vertically.  Other  variables  included  the  use  of  a  grain-nucleating  agent  in 
the  face-coat  of  the  shell  mold  and  superheat.  Some  of  the  molds  made  in  the  LaPorte 
foundry  were  bottom-filled  and  some  were  top-filled. 

Computer  simulations  at  The  University  of  Arizona  revealed  that  castings  with 
larger  grain  sizes  had  less  porosity.  Permeability  for  the  flow  of  shrinkage  metal  through 
the  mushy  zone  is  higher  for  larger  grains.  Hence,  pressure  within  the  mushy  zone  is 
greater  than  in  a  mushy  zone  in  a  casting  with  fine  grains.  The  greater  the  pressure  in  the 
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mushy  zone,  the  less  likely  porosity  develops.  Superheat  and  mold  preheat,  according  to 
the  simulations,  enhanced  progressive  solidification,  resulting  in  somewhat  less  porosity. 
Nitrogen  content  was  not  an  important  variable. 

The  horizontally  oriented  plates  cast  in  the  Whitehall  facility  had  extensive 
indications  of  porosity  as  picked  up  by  fluorescent  dye  penetrant  inspection  (FPI).  The 
“non-nucleated”  plates  had  less  porosity-indications  than  did  the  “nucleated”  ones,  and 
the  least  amount  of  porosity  was  found  at  the  intermediate  thickness.  The  porosity  in  the 
same  plates  was  also  analyzed  using  quantitative  metallography.  Plates  cast  in  molds  with 
a  nucleating  agent,  exhibited  porosity  of  less  than  0.05  vol.  %  to  about  0.2  vol.  %.  With 
the  use  of  a  nucleating  agent  in  the  mold-face  coat,  the  porosity  was  greater.  The  worst 
cases  had  0.2  to  0.45  vol.  %.  In  some  cases  the  porosity  was  not  uniformly  distributed, 
with  as  much  as  10  vol.  %  measured  near  the  surface.  A  weak  positive  correlation  ( = 
0.5  )  was  observed  between  quality  assessments  made  by  FPI  and  by  the  quantitative 
metallography. 

The  experimental  castings  made  at  the  LaPorte  foundry  also  exhibited  porosity, 
but  the  porosity  was  almost  entirely  internal.  Unlike  the  castings  made  in  Whitehall,  the 
porosity  was  not  surface-connected.  In  terms  of  incidental  near-surface  porosity  (not 
surface-connected),  castings  in  the  top-fdled  molds  displayed  significantly  higher 
amounts  than  the  castings  made  in  bottom-filled  molds.  Presumably,  with  bottom-filling 
there  is  significantly  less  turbulence  and  splashing  of  the  melt  as  it  fills  the  mold  cavity, 
which  results  in  less  nonmetallic  inclusions  that  nucleate  pores.  The  incidental  near¬ 
surface  porosity  was  associated  with  the  high  concentrations  of  internal  porosity, 
particularly  in  the  thinner  plates. 

In  addition  to  assessing  the  porosity,  grain  sizes  and  secondary  dendrite  arm 
spacings  (SDAS)  were  measured.  The  plates  cast  at  the  Whitehall  facility  had  average 
grain  sizes  that  ranged  from  0.1  mm  to  9  mm.  Castings  made  in  a  mold  with  a  nucleating 
face-coat  were  markedly  smaller  that  those  made  in  a  mold  without  the  nucleating  face- 
coat.  Even  with  the  nucleating  face-coat,  however,  some  of  the  thickest  plates  (0.120  in.) 
contained  relatively  large  grain  sizes  of  greater  than  1  mm.  The  average  SDAS  of  the 
castings  varied  between  22  and  37  pm,  depending  on  the  thickness  of  the  plate. 

The  set  of  castings  made  in  the  LaPorte  facility  comprised  castings  in  two  top- 
fdled  molds  and  two  bottom-fdled  molds;  the  former  had  an  inoculating  face-coat,  while 
the  latter  had  no  face-coat.  With  a  nucleating  face-coat,  the  average  grain  size  ranged 
from  approximately  1  mm,  near  the  ingate,  to  0.05  mm,  near  the  extremity  of  the  casting. 
The  average  grain  sizes  in  the  castings  made  without  the  nucleating  agent  in  the  mold, 
varied  between  0.75  and  1.2  mm  along  the  length  of  the  casting.  In  the  overall 
population,  the  finer  grains  (<  0.5  mm)  were  columnar  and  met  at  the  centerline.  The 
larger  grains  (>0.5  mm)  were  more  randomly  oriented  and  equiaxed.  The  average  SDAS 
in  the  castings  made  in  the  LaPorte  foundry  ranged  21.4  pm  to  25.1  pm.  Small  pores 
were  distributed  throughout  the  thickness  of  the  castings,  while  the  larger  and  more 
complex- shaped  pores  dominated  the  centerlines.  The  larger  and  more  complex  pores,  at 
the  centerlines,  formed  between  grains. 

Perhaps  the  most  significant  finding  was  the  extensive  surface-connected  porosity 
found  in  the  castings  made  in  Whitehall,  while  there  was  hardly  any  or  none  found  in  the 
castings  made  in  LaPorte.  This  suggests  that  the  period  between  when  a  mold  is  filled  and 
when  the  chamber  vacuum  is  broken  differed  in  the  two  facilities.  More  importantly,  a 
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systematic  study  of  the  formation  of  surface-connected  porosity  should  be  effected,  so 
that  controlling  processing  conditions  to  mitigate  its  formation  can  be  ascertained. 
Internal  porosity  is  tolerable  because  for  critical  applications  subsequent  hot-isostatic 
pressing  (i.e.,  hipping  )  is  part  of  the  production  cycle.  Hipping,  however,  does  not 
eliminate  surface-connected  porosity. 
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